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We find a unique direction in the CMB sky around which giant rings have an anomalous mean 
temperature profile. This direction is in very close alignment with the afore measured anomalously 
large bulk flow direction. Using Monte Carlo simulations, we estimate the significance of the giant 
rings at the 3cr level and the alignment with the bulk flow at 2.5cr. We argue that a cosmic defect 
seeded by a pre-inflationary particle could explain the giant rings, the large bulk flow and their 
alignment. 
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I. INTRODUCTION 

One of the key assumptions in modern cosmology is 
statistical isotropy. The detailed data from the Wilkin- 
son Microwave Anisotropy Probe (WMAP) provides an 
opportunity to test this assumption. Indeed many au- 
thors have studied this issue directly and indirectly using 
various approaches and claimed the existence of a num- 
ber of anomalies in the data (see e.g. [1-21] and [221, [23| 
for recent reviews). 

In this paper we propose another approach to test sta- 
tistical isotropy. We study how much giant rings in the 
Cosmic Microwave Background (CMB) sky deviate from 
random behavior and estimate the significance of the de- 
viation. In section II we define the rings score as a func- 
tion of the direction the giant rings surround and gen- 
erate a rings score map for the masked Internal Linear 
Combination (ILC) map. We show that the ILC rings 
score map has a clear peak. We estimate the giant rings 
in the ILC map as a 3<j deviation from ACDM. Moreover, 
we find that the giant rings are aligned with another re- 
ported ACDM anomaly [2J-[26| in the form of a large 
cosmic bulk flow. 

In section III we discuss a cosmological scenario that 
could explain the giant rings, the large bulk flow and 
their alignment. It is this cosmological scenario, which 
involves the effects of a pre-inflationary particle [23, l2S| •> 
that actually motivated us to look for these giant rings 
in the first place. Section IV is devoted to discussion. 



II. GIANT RINGS IN THE SKY 

We begin with the following question: Are there un- 
usual rings in the CMB sky? For reasons that will become 
clear shortly, we choose to focus on the largest possible 
rings, namely those that reside in a band of width P 
around = 7r/2 with respect to some direction specified 
by a unit vector n (see Fig. [1]). The band is symmetric 




FIG. 1: The score is calculated for rings in a band of width /3 
surrounding the equatorial defined by a direction h. 



and so the range of angles considered is 



(1) 



We denote the average temperature of an infinitesimal 
ring by T{9^ n) and the mean of the map by Tq and use 
the following rings score to detect unusual rings in the 
sky 



R{P,n) = 



d{cosO) f^{0,h), 



(2) 
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where T{0^n) = T{0^n) — Tq. This score is chosen since 
we are not looking to find any particular shape of T(6>, n). 
Rather, we are searching for the direction in which the 
rings deviate maximally from random gaussian behavior. 
For this we simply need to weigh correctly the contribu- 
tion of each infinitesimal ring to our rings score. This is 
the reason for the d{cosO) in the score. 

There are some issues one needs to deal with when 
working with actual CMB data. First, to have enough 
statistics in each infinitesimal ring, the rings cannot be 
taken to be arbitrarily small and the integral must be 
replaced by a discrete sum. In the results reported below 
for the 7- year ILC map (given in 1° resolution) we took 
dO -^ A6 = 3°, but have verified that the results are not 
sensitive to A^. Secondly, for obvious reasons we would 



like the results to be insensitive to Galactic foregrounds. 
Hence we use the KQ75 mask which removes 29% of the 
WMAP7 sky and calculate the quantity 



^di; 



^)= ^T2(i,n)M(i,n), 



(3) 



i=l 



where M{i^n) is the number of pixels in the z'th ring that 
survived the KQ75 mask cut and T{i^n) is the difference 
between the average temperature in the i'th ring around 
the direction n and the mean of the masked map. 

We would like to test the isotropy assumption of 
ACDM via the rings score. With this goal in mind j3 
cannot be taken to be too small since this will increase 
the chance that the direction favored by i?dis(/^,^) has 
no significant importance and is merely a statistical fiuke. 
However, due to the mask we are using we cannot take j3 
to be too large either. The reason is that as we increase j3 
the average size of an infinitesimal ring becomes smaller 
and so the ratio between the number of pixels we are 
masking and the pixels we are keeping when calculating 
Tii^n) becomes larger for generic values of n. 

Note that even if j3 is small R^isiP^ ^) does not approx- 
imate i?(/3, h) well when h points roughly perpendicular 
to the galactic plane, because the corresponding rings lie 
mainly inside the mask. To be on the safe side, we ignore 
directions for which more than 80% of the rings have less 
than 30% unmasked pixels in each ring. Thus using the 
cut sky eliminates a small portion of the map around the 
north and south poles. 

To verify that the defined rings score is not biased 
by the mask we used the HEALPix package to generate 
100,000 random maps and masked them with the KQ75 
map. We found that the peaks in the corresponding rings 
score maps are uniformly distributed on the sky and do 
not favor any particular area (see Fig. [2]). 

In Fig. [3] we plot the rings score calculated for the 7- 
year ILC map with the KQ75 mask for P = 7r/6, 7r/3 and 





7.61 [mKT 



FIG. 3: From top to bottom: The rings score calculated for 
/3 = 7r/6, 7r/3 and 7r/2, respectively. The score is calculated 
on the 7-yr WMAP ILC map (degraded to iVside = 64) with 
the KQ75 mask applied. Ignored directions are marked in 
black. 



7r/2. The dark regions near the poles are ignored for the 
reason explained above. We see that at around Galactic 
coordinates (/,6) = (276°, —1°) there is a distinct peak.-^ 
In fact, as illustrated in Fig. 21 even with the naked eye, 
the giant rings are visible. 

There are some interesting aspects to this rings score 
map in general and the peak at (276°, —1°) in particular. 
The location of the peak of R^isiP, ^) is fairly insensitive 
to the value of (3 (as long as it is not too small or too 



FIG. 2: A histogram showing the location of the maximum 
rings score (with /3 = 7r/3) for each of 100,000 randomly gen- 
erated maps using the KQ75 mask. The maximum scores for 
27% of the maps lie within the mask, which covers 29% of the 
sky. 



^ Up to effects of the asymmetric masking the rings score is sym- 
metric under a 180° inversion and so a symmetric partner peak 
appears as well. Taking an asymmetric range, however, the 
direction (276°, —1°) is more prominent. 




of astrophysical or systematic effects and not of cosmo- 
logical origin. Of course, we cannot completely rule out 
a similar explanation for the giant rings, but there is 
support for a cosmological explanation. First, the rings 
score maps for the V and W frequency bands are almost 
identical to that of the ILC (see Fig. [5]). Secondly, the 
ecliptic pole is located at (276°, —30°) and is ~ 30° away 
from both the main rings score peak and the secondary 
peak at (248°, -34°). 

Further support for the cosmological origin of the giant 
rings is the intriguing alignment between the direction of 
the rings and the direction of the large bulk flow reported 
in [25]. According to [25| the bulk flow on scales of about 
100 Mpc/h has a magnitude of \v\ = 416 ± 78km/s to- 
wards {l,b) = (282° ±11°, 6° ±6°). The chance of such a 
large bulk flow to happen in ACDM on such large scales 
is about 0.5% [25|. As plotted in Fig. [6l for /3 = 7r/3 
the rings score direction is (276° ± 9°, — 1° ± 7°) and the 
distance between it and the Bulk Flow direction is 9lg^ 
degrees. The probability of a 9° alignment between two 
random axes in the sky is 1.3%. 



FIG. 4: Top: The 7 year ILC map, smoothed to 3° reso- 
lution. The KQ75 mask is faintly superimposed on the map 
and the rings are marked around the dominant direction. Bot- 
tom: The same map rotated so that the dominant direction 
is placed at the center of the map. Incidentally, two of the 
Cold Fingers discussed recently by the WMAP team [23] (the 
ones that include what [23] refer to as Cold Spot I and Cold 
Spot IT) fall nicely inside the cold ring (55° < < 85°), with 
their hot counterparts falling inside two surrounding hot rings 
(25° < 6* < 55° and 85° < 6* < 115°). 



large). To be precise, when varying /3 from 30° to 165°, 
taking 10 evenly spaced values, the peak moves at most 
by 2°. Repeating this test on 10,000 randomly generated 
maps (masked with the KQ75 mask), only in 14 cases the 
peak was as stable. This implies a 3<j deviation from the 
statistically isotropic ACDM model. 

As is evident from Fig. [3] most of the signal comes from 
the largest rings. To quantify this we calculate for each 
map {/3 = 7r/6, 7r/3 and 7r/2) the score 



5(/3) 



Rr, 



R 



(4) 



where i^max is the maximum of the rings score map and 
R and a are the mean and standard deviation of the 
rings score map and get Sii^ciP) — 7.09,5.85 and 4.93, 
respectively. We estimate the significance by calculating 
S{I3) on random maps and find that for P = tt/G^tt/S 
and 7r/2 only 0.16%, 0.57% and 3.97% of the maps get 
a higher score, respectively. This implies that most of 
the signal comes from a narrow band around = 7r/2, 
and that as we increase /3 we increase the noise without 
increasing the signal. 

Several features breaking statistical isotropy have been 
found in WMAP data that turned out to be the result 




FIG. 5: i?dis(7r/3,n) for the foreground reduced 7-yr WMAP 
V (left) and W (right) frequency band maps (degraded to 
A^side = 64). The score maps are very similar to one another 
and to the ILC score map. In both maps, the peak is in the 
same location, (276°,— 1°) and its significance vs. random 
maps is ^ 1%. The location remains the same also when 
tested in higher resolution of A^side = 128. 



Overall, with current data we have two large scale 
quantities (one is large scale in terms of the CMB and 
the other in terms of Large Scale Structure) that are 
slightly anomalous that point roughly to the same direc- 
tion. In ACDM there is no correlation between the bulk 
flow and the rings score and so a fair point of view is to 
attribute the alignment between the two to a statistical 
fluke (which is not so rare - a 1.3% effect) and to argue 
that both features are not anomalous enough to challenge 
ACDM. 

In the next section we would like to offer a different 
point of view. We show that the scenario of [23, [28| nat- 
urally explains the anomalous bulk flow, the giant rings 
and their alignment. Taking the above results at face 
value, this scenario explains a one-in-a-million effect in 
ACDM. 




FIG. 6: The rings score direction (276° ±9°, -1°±7°) and the 
Bulk Flow direction (282° ± 11°, 6° ± 6°). The la distance 
between the two is 9lg^ degrees. Areas for la and 2a are 
shown. 



III. A POSSIBLE COSMOLOGICAL 
EXPLANATION 

The fact that the bulk flow appears to be too large 
is known for quite some time now. What is fairly new 
is that the shear and octupole moments associated with 
the bulk motion appear to be consistent with ACDM [255 
(see also 0). This suggests that these higher moments 
are generated by the standard ACDM power spectrum, 
while the overall bulk flow is generated by a non-ACDM 
faraway source. 

A model with exactly this feature (motivated by 31 1 
and earlier works [32|-[3J) was suggested recently in [23, 
|28|. There, we studied some of the cosmological imprints 
of pre- inflationary particles (PIP). We found that each 
PIP provides the seed for a giant structure (a spherically 
symmetric Cosmic defect - SSCD) whose gravitational 
potential is determined by the PIP in the following way 



*o(fc) = 



\H 



12^/7reA;3 



(5) 



k^a{t)H 



where e is the slow roll parameter, A = dm/d(j) — m^/e/2 
(with (j) being the inflaton field) and, as usual, the effect 
is evaluated at horizon crossing. 

For simplicity and concreteness it was assumed in [28| 
that dX/d(j) and Ug — I are negligible to find 



$(r,z = 0) = AClog(r), C = 1.09 x 10" 



(6) 



Relaxing these assumptions one typically finds 

<^(r)-r", (7) 

with |a| <C 1. 

Both dnj and ([7j) vary very slowly over large distances 
and therefore are quite different than typical potentials 
generated by the ACDM power spectrum. As a result, 
such a SSCD has distinct cosmological imprints [28]. In 
particular, it can induce a large bulk flow from far away 
towards its direction while having negligible effect on 
higher moments of the bulk motion. Hence it fits neatly 
with the observations of |25j . 

The CMB signal of a single SSCD (seeded by a PIP) 
is affected by its magnitude (A in the case of ([6])) and its 
distance from the observer, denoted by zq. Setting the 
magnitude to produce the measured bulk flow for each 
zo, we remain with zq as a single free parameter. The 
CMB signal is made up of competing contributions from 
the Sachs- Wolfe (SW) and the late integrated SW (ISW) 
effects and as was shown in [28], it should be detectable 
in the CMB in the sense that it is larger than the noise. 
However, near zq ^ 3 the two nearly cancel out, so that 
a SSCD located there can account for the measured bulk 
flow while adding a low, but detectable, signal to the 
CMB sky that would not immediately stand out as an 
obvious violation of statistical isotropy. 

So how does one search for the SSCD in the CMB data? 
In particular, we wish to find a way to tell apart the 
CMB signal of a SSCD from that of an unusually strong 
structure generated by the ACDM power spectrum. 

An overdense ACDM structure will induce a cold spot 
in the CMB sky if located at the last scattering surface 
and a hot spot if located nearby. Because of the unique 
large distance behavior of ([6]) (or (O), a SSCD will induce 
a more complex imprint that spreads all over the CMB 
sky (see [23]). The shape of this imprint is azimuthally 
symmetric and its profile depends on zq. Therefore its 
generic signature is a spot surrounded by rings. An in- 
teresting case happens where the contributions of the S W 
and ISW effects almost cancel out in the low multipoles 
(the cancellation happens at a different zq for each mul- 
tipole). This can lead to the disappearance of the spot 
for certain zq values. Since the potential falls slowly with 
distance, the fact that the circumference is maximal at 
= 7r/2 dominates and so a generic signature is the ap- 
pearance of anomalous rings around ~ 7r/2 from the 
location of the PIP. Focusing on these rings lowers the 
possibility that the score will confuse an atypical ACDM 
structure with a SSCD seeded by a PIP. 

Moreover, if indeed a single SSCD is responsible for 
most of the bulk velocity observed in [25| , then it should 
be located very near to the Galactic plane, where the 
small signal will be contaminated by Galactic fore- 
grounds much more than the large profile. Even though 
part of the signal predicted by the SSCD lies in small an- 
gles, limiting the search and focusing only on large angles 
yields a signal that is cleaner both with respect to Galac- 
tic foreground and ordinary ACDM effects. 



For these reasons we defined the azimuthahy invariant 
rings score the way we did in the previous section. It is 
designed to detect a SSCD seeded by a PIP. 

To verify that the score works as it should, we simu- 
lated random CMB maps that include the contribution 
from a SSCD located at some specific direction at a cer- 
tain distance, and checked if the direction with the maxi- 
mal rings score is indeed near the SSCD direction. This is 
done by calculating the SSCD temper atu re imprint from 
the SW and ISW effects (taken from [28]) and adding it 
to the randomly generated map (we calculate the effect in 
a multipole expansion and add the first non- vanishing 10 
multipoles a^os to those of the random maps). The map 
is then rotated so that the SSCD is hidden at a desired 
direction. Next, we calculate the rings score on the map 
and check whether this direction indeed dominates the 
map. In Fig. [7] we plot the percentage of random maps 
with a hidden SSCD where the maximal rings score is less 
than 8° away from the hiding location. We see that, as 
expected, it resembles the S/N graph in [28[ (Fig. 6(b)). 
A detailed study of the relation between the giant rings, 
other CMB anomalies and the PIP model will appear in 
a future work. 



IV. CONCLUSIONS AND DISCUSSION 

In this paper we reported on a novel unexpected fea- 
ture of the CMB sky - giant anomalous rings. The sig- 
nificance of these giant rings by themselves, much like 
other "anomalous" features of the CMB, is far from be- 
ing overwhelming as it is merely a 3cr effect. Moreover, 
much like the absence of large angular correlations [1, 2], 
our findings are weaker when considering the full ILC 
map with no mask - the peak in the rings score map re- 
mains aligned with the bulk fiow direction, but the score 
is no longer significant vs. random maps. This could 
either mean that this is due to the contamination of the 
unmasked data or that the feature is weaker than sug- 
gested by the masked map. With data from Planck we 
should be able to determine which possibility is right. In 
fact, since Planck is expected to clean much of the Galac- 
tic foreground, we should be able to include the small 



data as well and see if the signal increases. 

Estimating the significance of features in the data is 
tricky in general and in the case of statistical isotropy 
in particular. Indeed many of the reported large scale 
anomalies in the CMB that imply violation of statisti- 
cal isotropy were recently deemed as stemming from a- 
posteriori choices of estimators [29, 35-38] that amplified 
the significance of the results. Among the claims against 
these anomalies is that they surfaced from a search of 
oddities in the data with no independent experimental 
evidence or prior theoretical motivation. 

The giant rings are different in this regard. First, 
the search for them was motivated by a theoretical sce- 
nario which by construction violates statistical isotropy. 
Secondly, they are aligned with another large scale non 
CMB "anomaly" - the bulk flow. This increases, in our 
opinion, the chance that our flndings could eventually 
lead to a real challenge for statistical isotropy. For this 
to happen more data is needed. 

Luckily there are two clear predictions of our scenario 
that could be tested already with Planck data. First, 
the weak gravitational lensing of the CMB by the SSCD 
(assuming it has the magnitude required to produce the 
large bulk flow) is quite distinct [39] (again because of the 
long range gravitational potential it induces) and should 
be detected by Planck. Secondly, the measurement of 
peculiar velocities via the kinetic Sunyaev-Zel'dovich ef- 
fect should improve quite signiflcantly with Planck. This 
should enable testing the claims of [40|, which are based 
on data from WMAP, and determine if the bulk flow is 
indeed anomalous at even larger scales (which will in- 
crease the signiflcance dramatically), and in which di- 
rection it points. If our scenario is correct then as one 
increases the size of the survey, the usual ACDM effects 
should become more negligible compared to the SSCD 
effect and the measured bulk flow will be more aligned 
with the giant rings. 




FIG. 7: Percentage of detection (to within a distance of 8°) 
of a hidden SSGD at different redshifts zq by the rings score. 
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